Ultra-light Axions: Degeneracies with Massive Neutrinos and Forecasts for Future 

Cosmological Observations 
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A generic prediction of string theory is the existence of many axion fields. It has recently been 
argued that many of these fields should be light and, like the well known QCD axion, lead to ob- 
servable cosmological consequences. In this paper we study in detail the effect of the so-called string 
axiverse on large scale structure, focusing on the morphology and evolution of density perturbations, 
anisotropics in the cosmic microwave background and weak gravitational lensing of distant galax- 
ies. We quantify specific effects that will arise from the presence of the axionic fields and highlight 
possible degeneracies that may arise in the presence of massive neutrinos. We take particular care 
understanding the different physical effects and scales that come into play. We then forecast how 
the string axiverse may be constrained and show that with a combination of different observations, 
it should be possible to detect a fraction of ultralight axions to dark matter of a few percent. 



I. INTRODUCTION: THE AXIVERSE AND 
COSMOLOGY 



There is a widespread consensus that the Universe can 
be accurately described by General Relativity and the 
statistical physics of particles and fields. The quantita- 
tive model to arise from such a description can be used to 
accurately predict the gross features of the Universe (such 
as, for example, its expansion rate and the spectrum of 
the cosmic microwave background) and some detailed 
characteristics (such as, for example, the abundance of 
light elements, the number of relativist ic species, the den- 
sity of baryons and the curvature of space). In fact, over 
the past few years, cosmological observations have led to 
an ever increasing number of precision constraints on a 
variety of cosmological parameters that can describe the 
very early Universe, the evolution of the Universe over 
many orders of magnitude in scale and the current state 
of the Universe. Again and again, the favoured cosmo- 
logical model has passed observational tests with flying 
colours. It has been dubbed the concordance model [1]. 

The concordance model is incredibly robust yet, at the 
same time, it predicts a whole new sector of fundamental 
physics which has yet to be understood. In the concor- 
dance model, 96% of the Universe is dark, that is, doesn't 
interact or only interacts very weakly with electromag- 
netic radiation. This dark sector divides up into two 
components: roughly 20% of it is in the form of dark 
matter, a gravitationally interacting type of pressureless 
matter which can clump, and 76% in the form of dark 
energy, which is gravitationally repulsive and can drive 
the accelerated expansion of the Universe at late times. 
The overriding goal of modern cosmology is now to un- 
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derstand the dark sector using the tools of what has been 
called precision cosmology [2]. 

The focus in the quest for dark matter and dark en- 
ergy is to model and measure the large scale structure of 
the Universe. By this we mean the distribution of galax- 
ies as a function of redshift, the anisotropics of the cos- 
mic microwave background and the weak gravitational 
lensing of distant galaxies. The progress in measuring 
these different observables has been tremendous and will 
continue, with new experiments either being planned or 
coming online over this coming decade. 

Given the promise of observational cosmology, the 
hope is to not only understand the nature of dark en- 
ergy and dark matter but also explore other aspects of 
the concordance model and, in particular, measure in 
more detail other parts of the matter and energy content 
of the Universe. One obvious component of interest is 
the relativistic dark matter sector in the form of, for ex- 
ample, massive neutrinos. Another possibility is the exis- 
tence of ultra-light scalar fields that arise in some versions 
of string cosmology in what is known as the string axi- 
verse [3]. Both of these components are intimately tied 
to the dark sector and will be the focus of this paper. 



Additional Ingredients in the Hot or Relativistic 
Dark Matter Sector 



The natures of dark matter (DM) and dark energy 
(DE) are completely unknown, and theorists argue in- 
tensely even about their very existence, yet their rel- 
ative necessity in understanding precision cosmological 
data has made their gross properties in terms of ^a? 
and to a lesser extent the dark energy equation of state 
parameters wq and Wa^ become an integral part of the 
concordance model. In addition to these and the other 
established "vanilla" cosmological parameters, there is 
mounting evidence for additional cosmological ingredi- 
ents, possibly coming from multiple sources. 
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Observations of the cosmic microwave background 
(CMB) and large scale structure (LSS) in recent years 
have consistently found the need for excess relativistic 
enersv density ^- 

This excess radiation is param- 
eterised in terms of the effective number of relativistic 
neutrino species, A^eff,rei [^j-fllj. as: 
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where is the energy density in photons fixed by the 
CMB temeperature. Even within the standard model of 
particle physics, with no neutrino masses, A^eff,rei can be 
non-integer and greater than three if neutrino decoupling 
is non-instantaneous and the thermal neutrinos are par- 
tially reheated by electron positron annihilation [l2|. It 
is often stated that three massless standard model neu- 
trinos are best described by A/eff,rei ^ 3.04. Any increase 
from this is thought of as "extra" . 

The radiation density at big bang nucleosynthesis 
(BEN) is constrained by the light element abundances, 
but bounds vary depending on the treatment of astro- 
physical uncertainties, new physics scenarios and im- 
proved measurements. For example, [l3| has A/eff,rei = 

2.5^0 9, while [l4] allow up to AA^eff^rel = A^eff,rei-3.04 = 

1.39 at 95% credible interval. However, the BEN limits 
need not apply at the CMB or LSS scales, since late 
decaying particles may increase the neutrino abundance 
after BBN, but before or after CMB formation. Ac- 
cordingly, A^eff,rei IS taken as a free parameter in most 
cosmological parameter estimations, with best fit values 
from WMAP7 N.^^^ei = 4.34t^j^ [2], ACT N,^^,,, = 
4.6 ±0.8 [7], SDSS-DR7 A^eff,rei = 4.78t};?^ [5]. Moti- 
vation for such extra radiation density is lacking in the 
standard models of particle physics and cosmology, how- 
ever many theoretical extensions of these models provide 
clues. Such a situation is indeed expected in the "Freeze- 
in" mechanism of producing asymmetric dark matter, if 
the relic particles decay to neutrinos [l5|, [l6j . 

The observation of neutrino oscillations requires the 
introduction of neutrino masses, which constitutes a hot 
dark matter (HDM) component (for a review of mas- 
sive neutrinos in cosmology, see and for a historical 
review of HDM see [18]). However the nature of neu- 
trino masses and mixing remains unknown, and there- 
fore so also does the exact effect of neutrino masses on 
cosmology. The situation is further complicated due to 
the degeneracy, on certain scales and for certain observ- 
ables, between massive neutrinos and other cosmological 
parameters [19]. 

The situation with regards the measurement of mas- 
sive neutrino parameters using terrestrial experiments is 
summarised in some recent fits by Giunti [20] . These fits 
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seem to favour not only massive standard model neutri- 
nos, the absolute mass scale of which can only currently 
be determined by cosmology, but also the inclusion of one 
or two species of massive "sterile" neutrinos, which have 
a large mass splitting from the standard model neutrinos. 
The fits to this model and cosmology favour the standard 
model neutrinos being approximately massless, and the 
sterile neutrinos to have masses in the eV range. Sterile 
neutrinos are even more ambiguous cosmologically, since 
incomplete thermalisation allows the effective number of 
massive neutrinos, A/eff,mass7 to also take non-integer val- 
ues. Cosmological fits for sterile neutrinos are given in 
[14], and forecasts are made in |21|], while forecasts for 
the cosmological measurement of standard model neu- 
trino mass splittings are made in [22^. It is also well 
known that cosmology can potentially resolve the neu- 
trino mass hierarchy as being the "normal" heirarchy if 
the sum of neutrino masses, Sm^y < 95 meV, an accuracy 
within reach if many cosmological probes are combined 

a. 

The existence of neutrino masses, while in some 
views of the history of particle physics strictly "be- 
yond the standard model" (BSM), is at least well es- 
tablished enough, both experimentally and theoretically, 
that many would not class them as BSM at all. Their 
effects, as we have seen, are included in most cosmologi- 
cal analyses, too. Their known existence, combined with 
the fact that HDM cannot account for all the dark mat- 
ter, is definitive evidence that the dark sector is multi- 
component. Another theoretically well motivated, cos- 
mologically important ingredient that may be necessary 
to explain curious features within the standard model 
of particle physics, namely the strong CV problem, is the 
(QCD-)axion [2J-[39| (for reviews of axion cosmology, see 
[40|, |4l|). Depending upon their mass, axions can consti- 
tute the full range of dark matter "temperatures" , from 
cold through warm to hot: a true feast for Goldilocks. A 
model dependent coupling to photons can make them not 
really dark at all, and constraints can be derived on this 
from dimming of supernovae [42]. Hot, thermal axions 
can contribute to iVeff,rei/mass5 but their weak couplings 
make this contribution fractional [s^- For our purposes, 
the standard axion will be considered part of A/eff or 
appropriately. 

Axions in the mass range 0.7 eV < ^ 300 keV 
are excluded by cosmology for a variety of reasons 
[33, m, |43|. Sub eV mass axions contribute as HDM, 
and their mass is limited by constraints on A/eff: in exact 
analogy to neutrinos, they cannot be too heavy. Heavier 
axions, which are too heavy to be HDM, if they couple 
to photons, are restricted by their decays/inverse decays 
via effects on BBN, CMB distortions, and concordance 
between BBN and CMB determined values of the baryon 
to photon ratio. In addition, in this scenario, early axion 
decays to photons dilute the effective number of neutrino 
species, creating more tension with the large measured 
values of Aeff quoted above. 

Goldilocks properties of axions are abundant as they 
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make multiple changes in their dark matter temperature 
as one moves through their possible mass spectrum. One 
normally considers heavy weakly interacting massive par- 
ticles (WIMPs) with GeV masses as cold dark matter 
(CDM), and light neutrino-like particles with eV masses 
as HDM, and axions can indeed populate these masses 
and temperatures in the same way. However, very light 
axions with m < 1/ieV once again constitute CDM, and 
the types of limits given above from couplings to photons 
cease to apply [4^]. As we will see through the course of 
this paper, another transition occurs when these light 
CDM-like axions become lighter still and their quantum 
properties cause them to behave cosmologically like HDM 
again. This range of behaviour has to do with competi- 
tion between various physical processes which come in 
and out of dominance as coupling properties and relic 
density contributions of axions vary with mass and cos- 
mic evolution. We do not know the fundamental model 
and parameters that would exactly determine axion be- 
haviour. Names can also be deceptive and as we will see, 
the QCD axion is not the only axion relevant for cosmol- 
ogy- 

In addition, from BSM particle physics there are many, 
many, candidates for the CDM that differ from the stan- 
dard WIMP, the neutralino of the (minimally-) super- 
symmetric standard model, in lesser or greater ways. 
These differences may affect their properties during in- 
flation, during baryogenesis, or during BBN. These early 
time effects are of no concern to us here. As long as the 
dark matter is cold, we count it into Qc- Late time ef- 
fects, such as possible decay, or couplings within the dark 
sector, while cosmologically relevant, are also beyond the 
scope of this paper. 

From the point of view of non-linear structure forma- 
tion and cosmological phenomenology, dark matter self- 
interactions [44] and more novel ingredients like "Fuzzy 
Cold Dark Matter" (FCDM) [45] have been proposed as 
resolutions to the problems of cuspy dark matter halos, 
and the large predicted but unobserved numbers of dwarf 
galaxies in the standard CDM model (the well known 
"missing satellites" problem). This has lead some to gen- 
eral consideration of Bose-Einstein-Condensate (BEC) 
dark matter (see, for example, [4614481], and in the case 
of axions i50j). Indeed, the numerical simulation of 
[sij showed that the presence of such an ultra-light scalar 
condensate indeed reduces the number of dwarf galaxies, 
but in fact does very little to the cuspy density profile. 
However there are also many unaccounted for factors in 
standard galaxy formation models with CDM that may 
affect the formation of cusps and dwarf galaxies, such as 
baryon physics and supernova feedback (see, for example, 
(52I and references therein). 

A huge amount of research in modern cosmology goes 
into models for the 76% of DE that treat it other than as a 
cosmological constant i53L f or example modified gravity 
(5^ and quintessence j55l - l57( . Many of these models alter 
the equation of state, w{z)^ of DE, an effect degenerate 
with neutrino masses at some scales. In particular, mod- 



els where there is a component of dark energy with effects 
at high redshift (early dark energy, EDE) are known to 
share many degeneracies in their effects on cosmological 
observables with extra relativistic energy density, mas- 
sive neutrinos, and other forms of hot dark matter or 
any other structure-su ppre s sing cosmological ingredients 
(see, for example, [23|, l58|, l59|). The potential degen- 
eracies can, however, be broken by the use of multiple 
observables ^]. We will not be studying EDE in this 
work, but note that any potential detection of such an 
exotic component can only be truly qualified if all other 
aspects of cosmology with potentially similar effects are 
well understood. Finally, massive neutrinos are a key 
ingredient, along with a modified inflationary period, in 
allowing the model of [60] to fit the data and analyses of 
[61. §2] without the inclusion of a DE sector, so that the 
future success or failure of this non-standard model, too, 
must hinge on thorough understanding of the structure 
suppressing DM species. 



The Axiverse 

Cosmologists often invoke the existence of light scalar 
fields in the late universe as DM and DE components, for 
example in theories of quintessence, coupled quintessence 
(e.g. [63]), chameleons [64, 65], unified dark matter, and 
the Bose-Einstein condensates mentioned above. From a 
particle physics/string theory point of view these ingre- 
dients come up against two main problems: fifth- force 
constraints, and cosmologically light masses, which are 
not unrelated^. 

For a (coherent) scalar field to be cosmologically dis- 
tinct from CDM, or to play a quintessence like role, it 
must be very light: 10"^^ eV < < 10"^^ eV. Gauge 
invariance and Lorentz invariance then allow this scalar 
to multiply terms in the Lagrangian of the standard 
model fields, leading to problematic long range "fifth 
forces" . Unless a symmetry forbids them, these couplings 
should be universal and cannot be restricted ad hoc to 
the dark sector alone. If the field is to have an origin in 
new physics beyond the standard model then its light- 
ness and stability also become hard to explain without 
introducing additional hierarchy problems. For this rea- 
son, scalars in the late universe are considered generally 
problematic in models of particle physics; keeping them 
under control is one key motivation for moduli stabilisa- 
tion in string theory (see [69| and references therein). 

There is, however, at least one generic source of light 
scalars coming from high energy physics that evades all 
of the problems highlighted above: the so-called "String 
Axiverse" |3j] (where a more detailed version of the fol- 
lowing argument is given). Axions, as mentioned above. 



^ Addressing these issues for chameleons/string moduli has been 
looked at in e.g. [66l468[ . 
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were first motivated to solve tiie "strong CV problem" 
of QCD, where the problematic CV violating parameter 
occuring in the Lagrangian as £ D OF^iyF^^ is made 
dynamical as the Goldstone boson of a spontaneously 
broken global U{1) symmetry, and driven to its CV con- 
serving value by a potential induced non-perturbatively 
by QCD instantons. The quantum of excitation of is 
then the axion. The global symmetry is broken at the 
scale /a, whilst the non-perturbative physics giving the 
axion its potential switches on at a scale ja. This makes 
the axion a pseudo-Nambu- Goldstone boson (PNGB). In- 
deed, ultra- light axions and other PNGB's motivated in 
high energy physics, such as in [70., Jl], were studied as 
a solution to the non-zero cosmological constant prob- 
lem as early as 1995 [72], before the current vogue for 
quintessence in cosmology began with the supernova ob- 
servations of 1998 ^UM- 

We would like string theory to furnish us with the QCD 
axion and its solution to the strong CV problem. Ax- 
ions will always arise in string theory compactifications 
[75|,[765 as Kaluza-Klein zero modes of antisymmetric ten- 
sor (form) fields analogous to the Maxwell tensor, F^j^. 
These terms appear when the form fields are compacti- 
fied on closed cycles in the compact space, with 3-forms 
being compactified on closed 3-cycles, 2 forms on closed 
2-cycles etc. The number of axions arising due to the 
existence of a given form field is given by the number 
of closed cycles of the corresponding order; the relevant 
fields, however, are string theory dependent. All string 
theories contain the so-called "model independent" ax- 
ion arising from compactification of the antisymmetric 
partner of the metric, 5^^^, on closed 2-cycles. Generic 
string theory compactifications capable of realising re- 
alistic theories of high energy physics are highly com- 
plicated topologies, containing many hundreds of closed 
cycles (the source of the string landscape), and thus give 
rise to many axions. 

The underlying symmetries require axions to possess a 
shift symmetry, ^ -\- 27r^ and so their potential must 
be periodic. The Lagrangian for such an axion takes the 
following form: 



jC=i^{dOf-A^U{0) 



(2) 



where U{0) is some periodic potential. Bringing the ki- 
netic term into canonical form we define the field (j) = faO^ 
with Lagrangian: 



c = -id4>f-v{<p) 



(3) 



where V{(j)) is again a periodic potential. Expanding the 
potential in powers of 0//^, all the couplings of the field (j) 
come suppressed by the scale /o, and from the quadratic 
term we find that the mass is given by: 



A^ 

J a 



The symmetry breaking scale, fa and the scale of the 
potential. A, are both determined separately for each ax- 
ion, and depend on the action, due to non-perturbative 
physics on the corresponding cycle: 



f Mpi 



A 



(5) 



where Mpi is the reduced Planck mass: M^i = I/SttG. 
Solving the strong CV problem requires S > 200 [l, [76| , 
giving rise to stringy values of fa ^ 10^^ GeV, and this 
should be roughly constant for all these axions. The ex- 
act value of 5, however, scales with the area of the corre- 
sponding cycle (itself set by the scalar modulus partner 
of the axion), so that small variations in the area lead to 
exponential variations in the scale of the potential, and 
thus the axion mass. 

Hundreds of cycles of varying sizes therefore lead us to 
expect the appearance of at least some extremely light 
axions, given the following scenario. Axions generically 
get their masses lifted to high values at tree level, how- 
ever the required lightness of the QCD axion necessary 
to solve the strong CV problem, therefore avoiding such 
liftings of the mass, implies that other axions too may 
survive as light and stable. If string theory solves the 
strong CV problem by giving us the QCD axion, then 
the axiverse appears as a natural source of light scalars 
for cosmology. Some authors have been able to explicitly 
construct realisations of this scenario although 
the axiverse paradigm is expected to be much more gen- 
eral than these specific constructions. 

The axion shift symmetry also enforces that couplings 
to fermions appear derivatively as (90), leading to fac- 
tors of momentum at axion- fermion vertices. At low mo- 
mentum this suppresses long range forces on fermions by 
factors of k/fa and string axions avoid fifth force con- 
straints. In fact all axion couplings, including self cou- 
pHngs in the scalar potential, couplings to gauge fields 
like the photon due to higher dimensional operators, and 
topological couplings like the original OF^j^F^^^ come 
suppressed by this high scale. Therefore for cosmological 
purposes we will consider these axions to be completely 
decoupled, non-interacting massive scalar fields with po- 
tential: 
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(6) 



(4) 



Thus these axions are completely described by their mass, 
which we take to be a free parameter. We will con- 
sider dark matter axions with masses as low as m ~ 
10-^2 eV in the presence of additional CDM, neutrinos 
and a cosmological constant. Axion masses go as low 
as m < 10~^^eV, at which point the axion behaves as 
quintessence, as in the scenario of [72| described earlier. 
We will not consider quintessence axions. 

Some interesting features of the cosmology of string 
axions and their relation to infiation, the production of 
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gravitational waves and isocurvature perturbations were 
explored in [36|, [8Q| - [83| . String axions may also effect 
astrophysical phenomena, for example through black hole 
super-radiance (s^ . 

Ultra-light scalar fields are known to share qualitative 
features and many degeneracies in their effects on cos- 
mology with massive neutrinos and thus with many other 
cosmic ingredients [85] . It is the aim of this paper to ex- 
plore the effects of ultra-light string axions on the CMB 
and LSS, and thus these degeneracies, in detail. 



II. ULTRA-LIGHT AXIONS VS MASSIVE 
NEUTRINOS 

Light species of particles, such as massive neutrinos 
with rriu < 1 eV, can act as hot dark matter and suppress 
formation of large scale structure via free-streaming [86^ . 
Below the free streaming scale, kps^ which is determined 
by the temperature at which the species becomes non- 
relativistic, and therefore by the mass of the species, the 
HDM cannot cluster. During matter or A domination 
113: 



with these quantities again defined as in [90]: 



0.82 



(1 + ^)2 
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If there is a fraction of matter, /, in such a non-clustering 
species then the overdensities in matter grow as (5 ~ a*^, 

kps- This be- 
in the matter 



with q = l/4(-l + V25 - 24/) for k > 
haviour leads to the formation of "steps" 
power spectrum [85] . The size of these steps was first esti- 
mated in PI to be AP{k)/P{k) ^^8/^ (/^ = ft^/nm)- 
Fits for the steps can be found in [87] and [ss} . 

A qualitatively similar feature occurs in the presence 
of ultra-light scalar fields with ^ 10~^^eV, such as 
string axions, but the physics behind this process is quite 
different to the case of neutrinos or any other eV mass 
particles, such as the QCD axion d, EH, iH, [sij . 

A scalar field with the quadratic potential of Eq. [6l 
decomposed into homogeneous and inhomogeneous com- 
ponents as ^(^, r) = (I)o{t) + r) has the following 
equations of motion, to first order in cosmological per- 
turbation theory about a homogeneous, flat FLRW back- 
ground Hiii: 







bi + 2H^i + (m^a^ + k^ 



h = - ^4>oh 



(8) 
(9) 



where a is the scale factor of the FLRW metric, over- 
dots denote derivatives with respect to conformal time 
r, H = d/a and h is the scalar metric perturbation in 
conformal Newtonian gauge, as defined in [90j. The den- 
sity and pressure in the field are derived in the usual 
way, to first order, from the energy- momentum tensor. 
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(10) 

(11) 

(12) 

(13) 
(14) 



The equations of motion are oscillators, with oscilla- 
tions beginning when the mass overcomes the Hubble 
friction in Eq. [51 Using a WKB approximation we can 
solve Eqs. [HI [9] approximately and obtain the sound speed 
in the scalar field fluid perturbations, c,s = 6P/6p, as an 
average over the period of oscillation 



1; k > 2ma(i 



k < 2maa 



(15) 



On scales k > 2maa where the sound speed is 1 the per- 
turbations in the fluid are relativistic and overdensities 
will not grow, with the total overdensity in axions plus 
CDM scaling as (5 a*^, just like the case of massive 
neutrinos. On scales below this the sound speed goes to 
zero and overdensities grow as (5 ~ a, just like pure dust 
CDM. This behaviour is related to the Compton wave- 
length of the ultra-light particles, and implements the 
idea FCDM [H. 

There emerges a new scale, /Cm, analogous to the neu- 
trino free streaming scale, for ultra-light scalars. Modes 
with k > km enter the horizon whilst the sound speed is 
relativistic and will display a suppression of power, while 
modes with k < km enter the horizon as the sound speed 
is decaying to zero and cluster as ordinary CDM. These 
considerations give: 



1/3 



km ^ kp 



km 



1/4 



1/2 

— 1 ' k ^ k 



m, 
Ho 



(16) 



where Qm is the total fraction of the critical density in 
baryons, CDM, axions, and neutrinos (if they are non- 
relativistic at these scales), and Hq is the Hubble scale 
today. We note that this definition of km reproduces 
the same scaling with mass as the definition used in [^85^ 
where km is defined as the scalar field Jeans scale during 
matter domination evaluated at the redshift when scalar 
field oscillations begin. We further note that in our fits 
for the matter power spectrum made in subsequent sec- 
tions we will only be considering km < keq^ since massive 
neutrinos corresponding to the WMAP best fit values 
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have kps in this region. How much degeneracy there 
is between axions and neutrinos wih clearly depend on 
exactly how close km and kps are and how sensitive a 
particular observable in a particular survey is to physics 
on these scales. 

In the fits we use km described in [89| , where we add a 
bar to distinguish this fitted value, which fits the middle 
of the step, from the value derived above, which fits the 
start. The two differ by an order of magnitude for our 
fiducial cosmologies, which reflects the scale over which 
the transition in axion clustering behaviour occurs. The 
exponent 0^3 in these fits appears as ^m^ ^nd its fit- 
ted value of 0^3 = 0.4 is close to the expected value of 
as = 1/3. Finally we note that km enters the horizon 
when H ^ {H = aH)^ at exactly the same time when 
scalar field oscillations are expected to begin. This is 
again consistent with the definitions given in [85]. How- 
ever since Eq. [15] only holds once the fields have already 
begun oscillations we should take this as a warning that 
the expressions for km given by Eq. [16] will only be ap- 
proximate. 

Here we see the physical difference between suppres- 
sion of structure by ultra-light scalars and the free- 
streaming of neutrinos. Free- streaming is related to a 
change in temperature causing the particles to become 
non-relativistic when the temperature drops below the 
mass. Ultra-light scalars, if treated as particles, would 
still be relativistic today and their "free-streaming" scale 
would be larger than the horizon. However, ultra-light 
scalars form a condensate (see, for example, [49j and ref- 
erences therein), and we treat them as a classical field. 
As such it is the sound speed of perturbations in this 
condensate, which depends only on the mass and scale 
factor, not the temperature, which determines whether 
overdensities can form. This leads to an interesting co- 
incidence of scales: an ultra-light scalar with a mass in 
the range of 10~^^ eV will suppress structure on approxi- 
mately the same scale as a neutrino with a mass (9(10^^) 
times greater, which will make up about 1% of the total 
energy density. However, there is an extra parameter to 
consider for ultra-light scalars that comes from the pro- 
duction mechanism of their relic density and leads to an 
extra degree of freedom when considering their effects on 
the matter power spectrum. 

The relic density of massive neutrinos is fixed by the 
mass of each neutrino species, assuming standard model 
interactions. For one massive neutrino this is approxi- 
mately given by: 



93.14/^2 eV 



(17) 



where h is, and will from now on always be, the Hubble 
parameter defined as usual by Ho = 100/ikms~^ Mpc~^. 
We see that both kps and are fixed by the neutrino 
mass, rrij^. However, as mentioned in Section [H for the 
case of sterile neutrinos this is not necessarily true. Ster- 
ile neutrinos are thermalised with active neutrinos via the 
mixing process, but the degree of thermalization depends 



strongly on the masses and mixing parameters [92*] and 
therefore, just like the case of A^eff,reh ^eff,mass can be 
given non-integer values to parameterise this. For ther- 
malised standard neutrinos, since is fixed by the mass, 
we have kps = kpsi^u)- However for a sterile neutrino 
this is no longer true, since — l^i/(miy, A^eff,mass) 

Axions have two contributions to their relic density. 
First of all there is standard thermal production due to 
axion couplings to the standard model. However, just 
like the self interaction terms in the potential, all of these 
couplings appear in the Lagrangian suppressed by powers 
of ja- For large, stringy values of ja ^ 10^^ GeV these 
couplings are very small and the thermal relic density of 
ultra-light axions due to them is negligible. 

There is also a second, non-thermal production known 
as the vacuum realignment mechanism. The axion arises 
from the spontaneous breaking of the Peccei-Quinn sym- 
metry at the energy scale fa- At this scale the parent 
scalar field acquires a vacuum expectation value and the 
Goldstone boson, which is the axion, acquires a random 
initial value: the initial misalignment, 0^. Later, the 
field acquires its potential, V^((/>), due to non-perturbative 
physics and once the mass of this potential overcomes 
Hubble friction, H{zosc) ^ the field will roll and 
generate a relic density that depends on the initial mis- 
alignment. 

Until 

Zqsc the axion contributes negligibly to the en- 
ergy density as a cosmological constant, in contrast to 
a massive neutrino, which scales as radiation before be- 
coming non-relativistic. The effects on the redshift of 
equality by axions and neutrinos are thus not the same. 
Axions reduce the amount of matter only and do not af- 
fect the expansion rate while frozen; once rolling they 
have but a transitory effect on the background expan- 
sion away from ACDM as they go through their first few 
oscillations [89]. Massive neutrinos, or extra relativistic 
species, in our parameterisation reduce Qa compared to 
redshift zero, but also increase the amount of radiation 
at early times. This has a non- negligible effect on the 
expansion rate and leads to a markedly different effect in 
small scale CMB anisotropics between axions and neu- 
trinos, as we discuss in Section IHI El Fig. [1] shows the 
evolution of the energy density in various components. 

Using that Zeq = 2.5 x 10^ftmh^Q2.7^ where O2.7 = 
TcMB /"^'^ ^ and the neutrino sector consists of three stan- 
dard massless neutrinos, the relic density produced by 
this mechanism is approximately given by: 



= 8.4 X io-5/i-3/2ei. 7 (^-^^ 



— Q^m^^i '1 ^osc ^ ^eq 



1/2 



(18) 
(19) 



where (j) is dimensionlessly given in Planck units. 

For a quadratic potential, where we have explicitly 
broken the axion shift symmetry and are strictly work- 
ing with generalised ultra-light scalars, we are essentially 



7 




FIG. 1: Ultra-light axion and massive neutrino evolution in the background, nia = lO""^*^ eV, Qa = 0.1^^, ^eff,mass = 3, 
rriu = 0.1 eV. Left panel: evolution of the energy densities in massive neutrinos and axions compared to vanilla cosmological 
components. Right Panel: Contributions to the critical density, Qi = pi/3H^. 



free to choose to give us the desired relic density, be it 
large or small, for any axion mass. Therefore for axions 
the scale for suppression of power and the relic density 
are separately under control via the two parameters 
and in contrast to standard massive neutrinos where 
both are fixed by the mass, m^y. It is, however, pertinent 
to consider questions of fine tuning for this production 
mechanism. We direct the reader to the discussions of 
i82.,^,i8i,[9l. 

Another feature of vacuum realignment affecting the 
axion relic density comes from the periodicity of the ax- 
ion potential. The periodicity, if the shift symmetry re- 
mains unbroken, leads to the axion field having a maxi- 
mum value given by (j)max = ^^f^- Therefore it is clear 
that for axions below a certain mass it is impossible, bar- 
ring anharmonic effects in the potential, to produce 0{1) 
values for f^a, which leads to the existence of what the 
authors of [3| call "the anthropic window". However, in 
the spirit of cosmological parameterisation, this need not 
worry us. If our model requires a larger amount of ax- 
ion energy density at a certain mass than we can produce 
with < (l)rnax then we can view this in the same way as 
A^eff* it may be telling us that there are many species of 
axion in that mass range, where the masses cannot be re- 
solved. This may, however, require additional fine tuning 
within the axi verse. For the fiducial models considered 
for forecasts in this work, however, we do not saturate 
this bound, and so the shift symmetry is preserved. 

Standard model neutrinos come in three species, and 
each species should be massive, with some hierarchy and 
degeneracy structure between them. The cosmological 
detection of this degeneracy using weak lensing was dis- 
cussed in [22]. In the axiverse scenario we have multiple 
species of axion, with their own mass splittings. Naively, 
then, because of the qualitative similarity in their effects 
on cosmology, we may expect to account for any discrep- 
ancy between terrestrial measurements of Nj^^ rriiy and 
mass splittings with the values determined by cosmology 



via the introduction of ultra-light scalars. In contrast, the 
possible existence of sterile neutrinos and other relativis- 
tic relics may obscure the possible cosmological effects 
of axions and close this observational window on them. 
It is one of the principle aims of this paper to go some 
way towards addressing these potential degeneracies, and 
indeed we expect many of them to be broken by consid- 
ering multiple cosmological probes in the CMB and LSS, 
in the same way as degeneracies between neutrinos, dark 
energy and initial conditions can be broken in this way 
[iqI, [23I, [i^ • However, due to complications in forecasting 
for the effects of a varying axion mass and of the effect of 
a neutrino hierarchy splitting, we will leave the analysis 
of this particular degeneracy for a future work, and here 
focus in our forecasts purely on the density for a single 
species of axion, and on degenerate massive neutrinos. 

We summarise in Tabled the relevant scales of Zosci ^m-, 
and krn foi" the axions used in our fiducial cosmologies. 
We also quote keq{fax = 0), keq{fax) (where fax = 
and kpsi'^u^z = 0) for comparison. Note that although 
we always have km < keq^ this is not always the case for 
km- We also note that one's definition of Zqsc is somewhat 
ambiguous: does one define it from when = H or 
when rria = 3i^, or somewhere in between; when slow 
roll is broken, or when the oscillations have settled down 
to CDM behaviour? This leads to an 0{1) multiplicative 
factor of uncertainty. In particular, this makes Zqsc with 
rua = 10~^^eV potentially very close to Zrec ^ 1100. 



III. THE AXIVERSE AND COSMOLOGICAL 
OBSERVABLES 

As discussed in Section [iTl ultra-light axions give rise 
to steps in the matter power spectrum, P{k). Fig. [2] 
shows this effect on the cosmology of WMAP7 [H, with 
the introduction of a single axion species with fraction 
fax- Large axion fractions, disallowing variation of other 
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nia (eV) 


km 


(/iMpc" 






Zosc 


lQ-29 




0.0058 




0.0575 


350 


10-30 




0.0027 




0.0267 


74 


10-31 




0.0012 




0.0124 


15 


10-32 




0.0006 




0.0057 


2.4 


keq (^faa 


= 0) = 


0.0136/1 Mpc-^ 





keqifax = 0.01) = 0.0135/l MpC" ^ 

kFs{rau 0.055 eV, z = 0) 0.0451/i Mpc"^ 

TABLE I: Relevant scales for our fiducial cosmologies with 
fax =0.01. /cm is the scale at which structure suppression 
begins, given by Eq. [161 km is the location of the middle 
of the induced feature in P{k)^ fit for m loyi . Zosc is the 
redshift at which axion oscillations begin, which has an (P(l) 
multiplicative uncertainty. 




/c[/iMpc"^] 



FIG. 2: The matter power spectrum for three cosmolo- 
gies shown with current measurements from SDSS [95], and 
ACT [62]. We show first the WMAP7 cosmology (dashed 
black line). We also show two axion cosmologies, both with 
rria = 10"^^ eV: fax = 0.1 (solid blue line), and fax = 0.01 
(solid black line) , with all other parameters held fixed at their 
WMAP7 values. Both axion cosmologies have only a small 
effect on the CMB power spectrum, but are clearly distin- 
guished in their effect on P{k), with fax =0.1 clearly ruled 
out by the data. 



tion of the Boltzmann equations obtained from a mod- 
ified version of the publicly available code CAMB {qtI - 
[9^. Our modification introduces a module to deal with 
scalar fields having a quadratic potential with mass large 
compared to the Hubble rate, the bulk of which involves 
accurately fixing the initial conditions and background 
evolution in the presence of rapid oscillations, and inte- 
grating such oscillations accurately ^. 

We exactly numerically solve the evolution of the axion 
field, ^, the difficulty of which stops us exploring the re- 
gion of parameter space with rria ^ 10~^^ eV, suggested 
by [3] to be the most interesting region to look for unique 
step-like features in the power spectrum with a high pre- 
cision galaxy survey or 21cm tomography survey. We 
are also limited to studying a single axion field, however 
our results will show that in fact, since constraints from 
some observables are mass independent, this is not a lim- 
itation. Our technique makes no use of the approximate 
treatments of axion sound speed and averaging used in 
the analysis of [85]. In addition, the mass range that we 
study is the one found in [85] to have the most tightly 
constrained axion fraction, but also the range in which 
the approximations used are least sound. Future observa- 
tions will bound this regime even more tightly; making 
reliable predictions for high precision measurements in 
this important regime requires an exact treatment such 
as ours. 

Throughout this section we will use our physical intu- 
ition about the suppression of power caused by ultra-light 
axions, and the similarities and differences with respect 
to neutrino free-streaming, to try and understand our 
numerical results. Where possible, we will be guided by 
analytic fits, but stress that these are meant for qualita- 
tive purposes only, and have some limited applicability, 
which we discuss. Analytic fits are not used in our fore- 
casts. We emphasise that the figures and discussion of 
parameter variation in this section are meant only for il- 
lustrative purposes, and are not meant in any way as 
parameter estimation from existing data, nor do they 
necessarily reflect the fiducial models of our forecasts. 
A Markov- Chain Monte-Carlo analysis for parameter es- 
timation in this model will be the subject of a future 
work. 



parameters, can easily be ruled out at current sensitivity, 
while a small fraction of around 1% is indistinguishable 
from ACDM using the power spectrum alone, c.f. Fig. 1 
of [19]. We will see that the smaller fractions used in our 
fiducial models, while still indistinguishable from ACDM 
with a single observable, can be distinguished using pho- 
tometric redshift information and/or a combination of 
observables. 

In this section we discuss in detail the theoretical ef- 
fects of ultra-light axions on the various cosmological 
observables. The effects are explored both analytically, 
using the fits of i96i], and through numerical solu- 



A. The Matter Power Spectrum 

In [i^ we derived fits, Tax{k^ fax)^ for the shape of 
the steps in P{k) in a flat universe containing radiation, 
axions, CDM, and a cosmological constant. A, but no 



^ For more details, or a copy of the code, please contact us via 
e-mail. 
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baryons: 

fdTd{k^ fax) = fcTc{kj fax) + faxTa{kj fax) 

= fdTaxik, Z, fax)Tc{k, fax = 0) (20) 

Here, and throughout this paper, fi = Qi/Qd^ fi = 
fti/ftrn SO that fi = fifd- In [89], these quantities were 
equal since /^^ = 1 in the absence of baryons. The dif- 
ference between fax^ fax is important in the functional 
form of Tax^ as we will see below. We explicitly show the 
redshift dependence of Td arising from Tax^ which cor- 
responds to scale dependent growth. Whenever we drop 
redshift dependence, it is assumed that z = 0. 

The matter power spectrum is related to the transfer 
function by: P{k,z) = A{k)T^{k)Dl{z)^em Di{z) is 
the growing mode, given for example in [87], and A{k) is 
the primordial power. Therefore, the step in the matter 
power spectrum in the model of [89] was given by: 

In our numerical studies using CAMB the cosmology 
contains, in addition to CDM, axions and A considered in 
(89| . the other standard ingredients of baryons, and their 
coupling to photons, massless and massive neutrinos. In 
the presence of baryons we model the full matter transfer 
function according to [96] as: 

Tm{k) = fdTdik) ^ hn{k) (22) 

Td{k) is the total dark matter transfer function, including 
CDM, axions and massive neutrinos, if present, and Ti){k) 
is the baryon transfer function. The fitted baryon trans- 
fer function contains the gravitational effects of the cou- 
pling to dark matter through its dependence on the mat- 
ter fraction = Qd-\-^bi the sound horizon s, the drag 
epoch Zd^ the epoch of equality Zgg, the scale of equal- 
ity /Ceg, and the Silk damping scale ksuk- There is also 
a dependence on these scales incorporated into the dark 
matter transfer function. In the case of ultra-light axions 
that do not begin their oscillations until the matter dom- 
inated era, these scales should all be altered to account 
for the change in matter content during these epochs. For 
example, Zeq fc-^-bZeq- Since the gravitational effect of 
the dark matter has thus already been accounted for in 
the baryon transfer function, the step feature modelled 
by Tax{k^ z^ fax) should only multiply the DM transfer 
function, but with the weighting for axion effects coming 
in as fax = ^a/^m that is Td{k) Td{k,z,fax)' The 
total matter transfer function is thus given by: 

Tm{kj Zj fax) — fdTax{k^ Z^ fax)Tc{kj fax — ^)~\~ fh^hik^ fax) 

(23) 

where Tc(/c, fax) and T5(/c, fax) are fit by the formulae of 
[96^ , with relevant scales modified by the presence of the 
ultra-light component. The distinction between fax and 
fax is especially important in the explicit form of Tax 
from [H. 



Therefore the step in the power spectrum caused by 
an axion component in the presence of baryons is given 
by: 

where Tm{k^ z^ fax) is given by Eq. [23l Note that this fit 
is related to the fit of |88|] used to investigate the effects 
of massive neutrinos by AP{k)/ P{k) = T^^{k^ fax) — 1- 
In particular, the step size is defined by: 

^(/,,)= lim f^^ikja.) 

Given that the baryon transfer function goes to zero 
faster than the CDM transfer function as k goes to infin- 
ity, we have: 

S{fa.) = S{fa.) (25) 

where S is defined in [89| , which produces a smaller step 
than the case with no baryons, S{fax) > S{fax)' This is 
caused by axions making up a smaller fraction of the total 
matter than of the dark matter alone, i.e. fax < fax- 

The "naive" fit: fax{kjax) = Tax{k,fax), corresponds 
to axions suppressing growth on CDM and baryons 
evenly, with no account made for axion effects on the 
sound horizon, drag epoch etc. This, as expected, re- 
produces the small scale limit. The modified fit for 
Tax {k^ fax) incorporates changes to the sound horizon. 
Silk damping scale and drag epoch using the fits of [^965 
and thus qualitatively captures the deviations from the 
smooth fit due to distortions of the Baryon Acoustic Os- 
cillations (BAO), which are seen in the full numerical 
solution. However, the fits presented here end up overes- 
timating the total amount of power suppression by a few 
percent. 

A step in the power spectrum corresponds to a change 
in the ratio of small to large scale power, which can 
naively be mimicked by changes in other cosmological 
parameters, such as the tilt of the primordial power spec- 
trum. If measurements at large scales are poor, the effect 
can also be mimicked by adding more CDM, which shifts 
the power spectrum over to larger k. Isolating the unique 
effect of a structure suppressing species requires precise 
observations at the relevant scale, kps or km [19]. We are 
considering axion species varying in mass over many or- 
ders of magnitude, so have a correspondingly large vari- 
ation of the scale km- Power spectrum measurements 
have varying precision over this range of /c^, and so we 
expect different constraints on the axion fraction, and 
possibly different degeneracies with other cosmological 
parameters, for the axions of different masses, if km for 
the different species falls into regions of different accuracy 
in the survey. In particular, we should expect stronger 
constraints from galaxy redshift surveys alone on heavier 
axions with larger km- However, since all the axions we 
consider have km < keq^ where survey accuracy is at its 
lowest, this effect should not be significant. 
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B. Baryon Acoustic Oscillations 

The theory behind Baryon Acoustic Oscihations 
(BAO) and their effect on the matter power spectrum has 
been known since 1998 in the work of Eisenstein and Hu 
[qH , and they are an important cosm ologsj i cal t ool used in 
the distance ladder (see for example 1 100- 106], and most 
recently [107], and for a review see [108]). Measuring the 
BAO to high precision is thus a key goal in modern cos- 
mology Here we briefly discuss a method of extracting, 
and hence working definition of, BAO from the matter 
transfer function. 

We define the linear BAO as in e.g. 



..full 



[k) 



lin 



2^2 



m,no osc 



(fc) 



(26) 




where T^,fuii(^) is the matter transfer function for a 
certain cosmology, either numerical or analytical, and 
Tm.no osc(^) IS defined as the oscillation free, smooth 
transfer function taken as an n— node cubic spline of 
^m,fuii(^) at points /c^, i = l,...,n chosen empirically 
to get the best smooth fit. 

There is a small distortion of Tax{k, fax) iii a cosmology 
with baryons, away from its smooth form in the region of 
the BAO. We now aim to give some analytic understand- 
ing of the reason for this, and therefore predict how large 
we can expect any BAO distortions due to light axions 
to be. 

First of all, we return to our transfer function of Eq. 



Tm{kj fax) — fdTax{k^ ^, fax)Tc{kj fax — ^)~^fbTb{k^ fax) 

We see that while Tax is smooth and should thus be cap- 
tured by the spline in Bun, we have also introduced a 
dependence on axion fraction into the baryon transfer 
function, which is oscillatory. 

In Eqn. 21 of [96| it is clear that the amplitude of the 
oscillations in the baryon transfer function has a detailed 
dependence on all the cosmological scales: the Silk damp- 
ing scale, the scale of equality, the redshift of the drag 
epoch, and the sound horizon at the drag epoch, and all 
of these scales depend on the matter content at the time 
when they are relevant, and at all such scales the lightest 
axions were frozen and not contributing as matter. The 
dependence of the sound horizon on the redshifts of drag 
and equality further implies that varying these redshifts 
by including an exotic species will cause a variation in 
the period of the BAO, as we can see again from Eqn. 21 
of [96], or from the fit used in [100]. 

If we hold the total amount of dark mater fixed and 
introduce ultra-light axions with Zqsc < ^eq then all these 
scales are shifted relative to where they would be if the 
dark matter were pure CDM. The shift is simple to com- 
pute: we simply alter the matter content in the equations 
for calculating these scales by a factor (1 — fax) to ac- 
count for the frozen axions. The distortions in this case 
can be easily understood from [96] . We effectively change 



k [hUpc 



FIG. 3: The BAO (at z = 0) for two models: a fiducial ACDM 
model, and a model with axions having iria — 10~^° eV, 
fax = 0.1. We also show the BAO for the same axion model, 
but where the smooth transfer function is splined assuming 
ACDM. 



^rn{z) during specific epochs (l)o in the notation of [96|) 
while keeping the fraction Vti^/Vtrn (at z = 0) constant. In 
a cosmology with massive neutrinos we must further shift 
these quantities by increasing the density in relativistic 
species at the relevant epochs appropriately for the mass 
of the neutrino species: a more complicated effect. Note 
that such distortions in the case of massive neutrinos are 
not modelled in either [87] or [88], where baryon oscil- 
lations are not present: it is only the scale dependent 
growth due to massive neutrinos that is considered there 
(see Section mrC]) . 

In Fig. [3] we show the BAO calculated numerically 
including the effects of ultra-light axions with m = 
10~^^eV, such that km < keq^ keeping Qd constant, and 
fitting a smooth T^o osc(^) using a cubic spline. We see 
an overall suppression of BAO amplitude caused by the 
presence of axions, despite an increased ratio ^b/^d in 
the radiation dominated era. 

Some model dependence enters in our definition of 
^no osc{k). For our ACDM model, we find empirically 
the best k points for our cubic spline are k = 0.001 
and 0.029 < k < 0.369 with Ak = 0.05"^. The step 
in the transfer function caused by axions is a smooth 
feature, and is best captured by introducing two extra 
points into the spline at /c = 0.0081,0.02, where the 
first is our estimate for km in this cosmology. The fre- 
quency of the BAO is fixed by the choice of the set 
{0.029 < k < 0.369, AA: = 0.05}, since the spline an- 
chors Biin = 1 at these values, therefore we cannot see 
any change in the frequency in the figure. That the same 
set fits Tno osc{k) (by eye) both with and without axions 



This is slightly different from the values chosen in [101]. 
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FIG. 4: Change in the sound horizon at recombination as 
function of axion fraction as a percentage of its zero axion 
value, calculated using the formulae of [9y]. 



tells us that any change in frequency is small. To quan- 
tify this, we note that the frequ ency of the BAO is set 
by the sound horizon as sin(5A:) |lQQj . In Fig. |4]we plot 
the change in the sound horizon calculated using the for- 
mulae of [96] as a function of fax as a percentage of its 
zero axion value. Axions always increase the size of the 
sound horizon, but even for large fractions the change is 
only by a few percent. The increase is linear, a fact not 
at all obvious from the relevant formulae. 

We may ask whether the BAO can be biased by our 
choice of T^o osc{k)^ effectively the underlying cosmology 
that we assume. In Fig. [3] we show an example of such 
bias by choosing to instead fit a smooth transfer function 
using the k points necessary for a ACDM cosmology, with 
no extra points to fit the step. This leads to an increase 
in the apparent BAO amplitude on large scales in the 
first few oscillations. 

There is, according to our fits, a significant degeneracy 
between fax and Vtch^ . We can introduce more CDM 
to fix the epoch of equality while ultra-light axions are 
present. Fixing a flat universe, this will reduce Vt\ by 
some small amount, the effect on the power spectrum 
being only through the normalisation [96|, and invisible 
in ratios such as the BAO and Tax- Restoring equality 
also restores the Silk damping scale and the drag epoch 
to their axion- free values, and this restoring effect dwarfs 
any small changes to the BAO through the alteration 
of the matter fractions fd-, fh- In this case the BAO 
distortions can be removed, and the overall suppression 
of power due to axions. Taxi also reduced as the in- 
crease in CDM shifts the power spectrum over to larger 
k. Avoiding such an alteration to the BAO will, however, 
have obvious effects on the CMB through changing the 
entropy per baryon, and will need further compensation 
for example in fitting the Hubble expansion rate by low- 
ering the DE equation of state, w. We discuss some of 



these issues in Sections IIIIDI and IIII E[ 

Real BAO measurements depend on redshift, can mea- 
sure the expansion rate as a function of which can 
break some of the simple degeneracies discussed here. 
In addition, their measurement in models including ax- 
ions will therefore be complicated by the scale dependent 
growth introduced by the z dependence of Taxi which we 
now discuss. 



C. The Growth Rate 

The growth rate is another useful cosmological observ- 
able, which has been measured by [109-112], and is used 
particularly in studies of modified gravity. It is defined 
as: 

a ma tio 

In standard ACDM it is known to be approximately scale 
invariant and behave with redshift as / = Vtrni^V^ (this 
is a useful approximation, but see for example [113] for 
a recent discussion). We expect the growth rate for cos- 
mologies including axions or massive neutrinos to pick up 
some additional scale dependence related to the appear- 
ance of steps in the matter power spectrum. 

In Fig. [5] we plot the growth rate as function of z 
for a standard ACDM cosmology, and for a cosmology 
with a fraction of axions fax = 0.01. The growth rate 
with axions is plotted at two different /c- values: the high- 
est and lowest from P(k) being kmin ^ 10~^/iMpc~^, 
kmax ^ 0.3/iMpc~^. The true growth rate with axions 
(insert of Fig. [5]) is seen to contain rapid oscillations that 
obscure the details of any step in f{k) at a given redshift. 
These oscillations are not observable, and should be av- 
eraged over to be consistent with the interpretation of 
scalar fields behaving as dark matter in the background, 
which is also only true on average as seen from the WKB 
solution to Eq. [51 To account for this easily we present 
an approximation scheme that is valid for both axion 
and massive neutrino effects on the growth rate when 
the fraction in these species is small, as is true for all of 
our fiducial models. 

We split the density into two pieces: pc+6 in the dom- 
inant CDM and baryon components, and pa+v in the 
sub-dominant axion and massive neutrino components, 
and do the same for the perturbations. Using that Pa+z^ 
is of order a few percent of we expand the expres- 
sion for / and drop terms of order p^+zy and pa+v^P as 
being second order. Next we note that the over densities 
will also be predominantly made up of CDM and baryons 
and so further expand in powers of 5pa+v / ^Pc+h and take 
terms of order {5pa+v / ^Pc+hY i pa+v{5pa+v / ^Pc+h) and 
(5pa+zy/^Pc+6 second order. The resulting expression 
for the growth rate is given by: 
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FIG. 5: The (reduced) growth rate (Eq.[28]), f{z), for rua = 
1Q~ ^° eV, fax = 0.01, with current measurements of [l09- 
[TTil- Red: f{c+h){k - kmax); blue: f(^c+b){k - kmin); black: 
ACDM. Insert: zoom in to an arbitrary region of z showing 
oscillations in the true growth rate (dashed lines) about the 
reduced growth rate. The reduced growth rate is seen to be 
a good qualitative tracer of the average, following the shape 
in z and demonstrating scale dependence in its amplitude. 



where f{c+h) is the growth rate in only the CDM and 
baryon components, but calculated in the cosmology in- 
cluding the exotics. We will call this the reduced growth 
rate. The true growth rate is a small perturbation about 
this, with all potentially oscillatory contributions iso- 
lated. The effects of the axions and neutrinos contribute 
to this reduced growth rate only through the gravita- 
tional couplings via the potential /i, and through the 
background expansion coming from 1-L. 

It is this reduced growth rate that was shown in Fig. [5l 
In the blow up insert, we see that as expected it traces 
somewhat the average of the oscillations in the true 
growth rate, and shows scale dependence. 

To more clearly show the scale dependence, in Fig. [6] 
we plot the reduced growth rate, /(c+6) as a function 
of /c, normalised to one on the largest scales so that 
z-dependence and normalisation to ACDM are absent . 
There is a clear step occurring at k i~ the same scale 
as the suppression of power in P{k). 

We have checked that the same growth rate is recov- 
ered in both the case where massive neutrinos are in- 
cluded exactly^, and where an averaging is taken over 
the oscillations present with an axion component. The 
accuracy of the approximation is to within a few percent 
of the total step size, which is itself only a few percent of 
the total value of the growth rate in a pure ACDM cos- 
mology. It is thus this reduced growth rate that we will 
use in our forecasts for galaxy redshift surveys (GRS), 
which require the growth as a function of /c in a redshift 
bin, and is output for them from our modified version of 



^ In this case the growth rate contains no visible-by-eye oscilla- 
tions, and is computed exactly using a numerical derivative of 
5. 



FIG. 6: The reduced growth rate, /(c+b)(/c), for three fiducial 
cosmologies: ACDM; ACDM with axions: ma = 10-^^ eV, 
fax = 0.01; ACDM with massive neutrinos, m^ = 0.055 eV, 
A^eff,mass = 3.04. Wc uormallsc all growth rates to one on the 
largest scales to account for normalisation by ACDM on the 
largest scales at arbitrary redshift. The results shown in solid 
are numerical, and the dashed line is that predicted by our fit 
for A/ in the presence of axions, Eq. 1291 The vertical dashed 
line on the plot shows the expected value of A^m, fit using the 
formulae of [sol], while the horizontal dashed line shows the 
expectation of Eq. 1301 

CAMB. In these forecasts we therefore make no use of 
the fitting / = l^^(z)'^^. 

Now we turn our attention to understanding the size, 
shape and position of this step in the growth rate in more 
detail. We have already described the effect of axions on 
the matter power spectrum using the transfer function 
Tax- We use the same transfer function for the overden- 
sities, i.e. 5m = Tax{k)S\ci:>M' Note that this is smooth 
as a function of z: Tax is already an on- average fitting. 
Note also that this is an alternative to parameterisations 
of scale dependent growth used in discussing massive neu- 
trinos in [23|, l87|. Substitution into Eqn. [271 immediately 
yields: 

1 f 

r r . -'- ax 

J - /(ACDM) -^TJTf— 
^ J-ax 

= /(ACDM) - ^f{k,z) 

(29) 

Using our fitting formulae we can calculate A/. We 
work in the regime where fax is fixed as a function of 
time, i.e. the axions have completed their transition to 
matter like behaviour, which will always be true for our 
fiducial models in the redshifts of interest, and in any case 
the fitting formulae break down where this condition is 
not satisfied. 

In the redshifts of interest, the z dependence of the 
additional term is mild relative to /(acdm)^ and so the 
shape of the growth rate as a function of modulo the 
oscillations, is largely unaltered by the presence of axions. 
In Fig. Owe have already shown the reduced growth rate 
as a function of and see that this is indeed the case. In 
Fig. [6] we also show the fit of Eqn.[29l Notably, since the 
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fit gives us tiie total growtli, /, ratiier tiian tiie modified 
growtii, tiie fit overestimates ttie size of tiie step. Tiiis 
can be taken as an indication of tiie amount by wiiicii 
tiie modified growtii is periiaps an underestimate of tiie 
true step. 

We estimate tiie ciiange in growtii rate amplitude at 
scales below km to be: 



0) = {l-q)^ -fa 



(30) 



for small fax- 

The most important feature of tiiis analysis lias been 
the identification of smooth, step-like, scale dependent 
growth occurring in models containing an ultra-light 
scalar field caused by the same physics as causes the sup- 
pression of power in the matter power spectrum, and also 
mimicking the corresponding effect due to massive neu- 
trinos, though to a lesser degree, as seen in Fig. [6l While 
the size of this effect is not accurately estimated ana- 
lytically, it's location in /c-space can be predicted with 
reasonable accuracy. The location of the step depends 
only on the mass of axion, while the size depends only 
on the density fraction. Measurements of the growth rate 
amplitude will therefore measure the fraction in axions, 
but be insensitive to the mass unless scale dependent 
growth can be resolved at km- Suppression of the growth 
rate, if measured in a non-scale dependent way, is clearly 
degenerate with the total total matter content, evident 
from the simple fitting in ACDM of / = rtm{z)^^ ^ and 
suggesting further positive correlation between fax and 
Qch^. Again, using a single observable, precise scale de- 
pendent measurements are required to isolate a unique 
signal. 

Weak lensing tomography also measures the growth 
rate, and we can use this to put tight constraints on the 
existence of smooth components, such as axions at scales 
smaller than km, or massive neutrinos at scales smaller 
than kps, from the amplitude change relative to ACDM 
jll4']. In the case of weak lensing, which measures the 
growth rate more accurately and via different means, no 
approximation is made at all and the full numerical evo- 
lution of the overdensity is used. In this case, the oscilla- 
tions are naturally smoothed by the redshift integral (see 
below). 



D. Galaxy Weak Lensing 

Galaxy weak lensing was first observed in 2000 by |ll5l - 

Ill8| . The measurement of the weak lensing power spec- 
trum is a direct probe of the dark matter. Through 
tomography, made possible by measuring galaxy pho- 
tometric redshift s, we also gain information in the ra- 
dial or temporal direction, and probe the growth rate 
and distance-redshift relation. This measurement of the 
growth rate gives constraining power for the presence of 
ultra-light axions. 



We first review some weak lensing basics. The ef- 
fect of weak gravitational lensing is usually split into 
two components: the complex shear (71,72) and the 
convergence with the shear being derivable from the 
convergence (see ^119j for a complete review). Let us 
denote x the comoving distance and r(x) the coordi- 
nate distance, defined by r(x) = sin (i^^/^x) fo^" 
a closed universe, r(x) = X foi" a flat universe, and 
^(x) = (—i^T) "^/^ sinh ((—i^r)-^/^x) foi" an open universe, 
where K is the curvature. 

The convergence in a given direction ft of th e sky is 
given by an integral along the line-of-sight |l20| 



(31) 



K{ri,x)= / W{x')S{x')dx' 
Jo 

where 6 is the density perturbation and 

W{x) = ln^H^g{x){l + z). 
is a weighting function, ^(x) is given by 



gix) = r{x) / 

Y 



dx'^(xO 



where n(x) is the source distribution: 



n{z) 



3 f z 



2^0 V^o 



(32) 



(33) 



(34) 



We can now expand the convergence in multipoles hiim, 
and define the convergence power spectrum by the rela- 
tion 



{f^lmf^l'm') — ^mm' Pi -> 



(35) 



and using equation (|3T]) under Limber's approximation, 
we have the following expression: 



Jo 



dx-^P(VKx),^) 



^^(x) 



(36) 



where Xoo stands for x(^ and P is the matter 

power spectrum. 

In this way, the convergence power spectrum depends 
directly on the matter power spectrum at redshift and 
thus tomographic information will measure the growth 
rate. The scale at which the growth rate is measured 
depends on the pixel size, and for our purposes this will 
measure below either km or kps for all axion and neu- 
trino masses considered. Therefore, the constraint on 
these species will be independent of this scale and will 
depend only on the density fraction, i.e. for the case 
of ultra-light axions where the density fraction depends 
only on the initial misalignment angle constraints on fax 
from weak lensing are expected to be independent of ax- 
ion mass^ but due to their similar effects will be degen- 
erate with the neutrino density fraction, and the dark 
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energy equation of state, the stronges t con straint being 
from the amphtude of the growth rate |ll4j . In addition, 
our method for calculating the weak lensing tomography 
comes directly from integrating the density parameters 
inside CAMB, and so is not limited by the approximation 
of the reduced growth rate used for GRS, as mentioned 
already above. 

Note that this effect in the growth rate is distin- 
guishable from many DE effects, which occur predom- 
inantly through modifying the distance-redshift relation 
via changes in the expansion rate: 

X{z)= r dz'H{z') (37) 

^0 

As we have already seen, axion and neutrino effects on 
the expansion rate during the matter era are small to 
non-existent. Establishing the distance-redshift relation 
in addition to the growth rate using lensing tomogra- 
phy serves to break degeneracies occurring between ax- 
ion/neutrino components and DE. For example, as dis- 
cussed earlier one can change the DE contribution to re- 
store keq and the canonical shape of P{k > keq), but this 
then demands a change in w if we are also to restore the 
Hubble expansion, leading to an apparent degeneracy be- 
tween w and m^y or fax^ which can be broken by lensing 
tomography or BAO measurement s_ that pin down the 
expansion as a function of redshift [94]. 

Because of the dependence on P(/c,2:), all of our in- 
tuition of the similarity of axion and neutrino effects 
on the matter power spectrum should carry over into 
galaxy weak lensing. Only our Fisher matrix analysis 
will show what degeneracies really exist in the fiducial 
models we investigate under the precision of the obser- 
vations in question. 

In Fig. Owe show the effect in the convergence power 
spectrum of adding an axion fraction fax = 0.01 in a 
species of mass = 10~^^ eV in various different red- 
shift bins. The unclustered species causes a suppression 
of power on small scales, which increases with increasing 
fraction in that species, as one would expect from the 
structure suppression in P{k). The errors on the con- 
vergence power spectrum expected form a large future 
survey are also shown (see Appendix [A] for survey param- 
eters). With a single redshift bin, these are not strong 
enough to distinguish a ACDM model from a model with 
axions at la. The amplitude of suppression of power is 
independent of axion mass, and like effects in the growth 
rate is smaller than the same effect due to massive neu- 
trinos. 

A possible uncertainty on constraints from weak lens- 
ing lies in the treatment of the non-linear regime, neces- 
sary for weak lensing calculations. Inside CAMB we use a 
standard version of halofit [121], which is optimised from 
simulations to standard CDM and may not accurately 
describe the non-linear effects of axions or massive neu- 
trinos. However, mistreatment of non-linear effects likely 
leads to an underestimate of the suppression of power on 
small scales [23], and so our projected weak lensing con- 
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FIG. 7: The convergence power spectrum for three redshift 
bins, with expected la errors from a large future weak lensing 
survey, defined in Appendix El Overlapping regions of error 
are shown in grey. The models compared are ACDM, and 
ACDM+/a,, with rua = IQ-^^eV, fa. = 0.01. 



straints should be conservative. Hannestad et al [94] es- 
timate that the effect of assuming no non-linear effects in 
neutrinos leads to an uncertainty in lensing observables 
of around 0.1%. For a discussion of the non-linear effects 
of neutrinos see |l22| , and for an N-body simulation with 
ultra-light scalars, see [5l|. 



E. The CMB 

1. Temperature Power Spectrum 

Here we discuss axion effects on the CMB temperature- 
temperature auto correlation power spectrum, Cj^ ^ the 
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TT pow er spectrum (for a review of CMB anisotropies, 
see |l23| , and for sp ecific applications to HDM and neu- 
trinos see [T23 . fT25| ). 

In Fig. [8] we show the effects of various axion fractions 
and masses on Cj^ . Axions have the effect of shifting 
the acoustic peaks by changing the epoch of equality. 
The changes in equality in these models, Azeq{fax) = 
Zeq{0) - Zeqifax). are Azeq{0.1) = 328 and A^eg(O.Ol) = 
32.2 respectively. 

The most pronounced effect for heavier axions occurs 
in the Integrated Sachs- Wolfe (ISW) region at very low 

as pointed out in [85]. As also pointed out in (sBj . 
this ISW effect is maximal for masses at the larger end 
of those we study, suggesting potentially strong CMB 
constraints caused by this. These heavier axions, as dis- 
cussed earlier, are making their transition from cosmo- 
logical constant to CDM behaviour very close to CMB 
formation. As such, they are contributing a significant 
fraction to the energy density while in a non-standard 
stage of evolution, and thus causing large variations in 
the potential on the largest scales where this transition 
is beginning. At high i there is no difference in the effect 
of axions of different mass, so small scale CMB measure- 
ments will constrain the total fraction possible in ultra- 
light axions independent of mass. Heavier axions than 
those we consider here that have completed their transi- 
tion to CDM in the radiation dominated era should be 
expected to have no effect at all on the CMB, as they are 
gravitationally indistinguishable from CDM. 

In addition, the low--^ ISW effect distinguishes axions 
from massive neutrinos. It is in the high-£ region of the 
TT power spectrum that evidence from the CMB for ex- 
tra relativistic energy density comes, due to the effect 
of this energy density on the pre-recombination expan- 
sion rate and corresponding increase in the amount of 
Silk damping leading to the required lower fluctuation 
amplitude [3, Il26| . The best way to explore the effect 
of extra relativistic species on the TT power spectrum 
is to hold fixed a number of quantities and isolate the 
extra effects at high-^. We can apply the same logic to 
further understand the effect of ultra-light axions since 
they affect the epoch of equality in the same way as ex- 
tra relativistic energy density: increasing A/eff,rei means 
more radiation, making equality later, while introducing 
axions with Zqsc < ^eq reduces the amount of dark mat- 
ter at equality, giving the same effect [125]. The amount 
of baryons, p^, should be fixed to keep the even and odd 
acoustic peaks in the correct ratio of amplitude. The 
epoch of equality, Zgq, being tightly constrained by the 
position of the first peak should then be fixed by changing 
pd. The angular size of the sound horizon, Os = Ts/Da^ 
should be fixed by changing Hq. Finally, the height of 
the first peak, which isolates the low and high-^ effects, 
should be fixed by changing the scalar amplitude, Ag. 

When all of these changes are made, the result is to 
completely remove the effect of axions from the CMB, 
while the same changes made in a cosmology A A/'g/ / > 
leave the expected effect in the Silk damping tail This 



is because the extra relativistic energy density has left 
an imprint in altering the expansion rate during the ra- 
diation era. In an axion cosmology, restoring equality, 
as mentioned above in our discussion of the BAO, also 
restores the Silk damping scale. Axions in the CMB are 
degenerate with a combination of changing pd and Hq, 
while neutrinos are not. These observations are impor- 
tant, since they will apply also to massive neutrinos: the 
CMB breaks the degeneracy between axions and neutri- 
nos due to neutrino effects on the expansion rate in the 
radiation era, and axion ISW effects in the matter era. 

E-mode CMB polarisation auto-correlation power 
spectra , Cf^, and temperature-polarisation cross- 
correlations, Cj^^ the EE and TE power spectra, are 
also useful cosmological probes, and they too are used in 
our forecasts. Polarisation will also be particularly im- 
portant if ultra-light axions have the model dependent 
coupling to E ■ B. In this case the polarisation angle 
can be rotated by A/3 ~ 10~^, an effect within reach 
of Planck and CMBPol. If an observation of this ro- 
tation were combined with other axion observations as 
described above, this would provide very strong evidence 
for the existence of axions with this coupling, but since 
the coupling is model dependent, polarisation rotation 
need not accompany other axion effects. Observation 
of polarisation effects without other axion effects would, 
however, rule out axions as a cause of the rotation. We 
do not include the possibility for such a rotation in our 
models. 

The CMB is an important constraint on another effect 
of axions, which we do not consider in this work: axion 
isocurvature perturbations. Axion isocurvature contribu- 
tions are significant if the energy scale of inflation is high. 
The authors of [3] assert that a low scale of inflation is 
natural in the axi verse scenario. If axions are not to over- 
produce isocurvature perturbations, the energy scale of 
inflation must be so low that gravity waves would not be 
significantly produced. The authors of [36^ propose this 
as a direct method of falsifying the existence of string 
axions, while the authors of [82, 83] propose using infla- 
tionary parameters in conjunction with axion isocurva- 
ture contributions as a method of quantifying fine tuning 
within string axion models. We propose to investigate 
axion isocurvature in detail in a forthcoming paper. 



2. CMB Lensing 

CMB lensing has recently been detected and measured 
by ACT [61, 127], and its observation by Planck will be 
a powerful cosmological probe. 

CMB lensing, like galaxy weak lensing, is a direct 
probe of the DM density and expansion history between 
the surface of last scattering and us. Hence, in terms of 
CMB observables, the EDE fraction and Sm^y are most 
strongly constrained by CMB lensing [23]. However for 
our purposes, we cannot include CMB lensing and galaxy 
weak lensing together since addressing correlations be- 
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FIG. 8: Effect of axions on the CMB power spectrum. Left panel: TT power spectrum. Right Panel: magnified low I region 
showing the ISW effect. 
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FIG. 9: The CMB lensing power spectrum. 



model our surveys to be similar to the proposed Euclid 
mission. We combine these results with Fisher forecasts 
for an all sky CMB survey based on the Planck mission. 
The Fisher matrix provides the lowest possible intrinsic 
statistical uncertainty, thus results from Fisher forecasts 
tend to be optimistic when compared to real results. We 
describe in detail our own implementation in galaxy red- 
shift surveys and galaxy weak lensing tomography, while 
our CMB forecasts are made using the pre-existing pack- 
age "FisherCodes" . 

Fisher matrices are widely used in cosmology to fore- 
cast uncertainties in parameters O of an underlying 
model of an observable x given an uncertainty in a mea- 
surement of X, (Jx- If we assume that the uncertainties in 
our model are Gaussian, the Fisher matrix can be calcu- 
lated directly from derivatives of the model, and is 



tween these observables is beyond the scope of this pa- 
per. We do not use CMB lensing in our forecasts, but we 
show the power spectrum here and make some comments 
for the sake of completeness. 

We show the lensing power spectrum in Fig. [9l A large 
fraction in axions produces a noticeable effect in reducing 
the lensing power. However, if as above we fix Zeq and 
Os again the effect can be totally removed. This occurs 
due to the degeneracy of axions with other parameters 
in their effects on the matter power spectrum, as already 
discussed. 



IV. COSMOLOGICAL OBSERVABLES: 
FORECASTS 

In this section we consider the potential to identify 
traces of axion dark matter with some of the cosmolog- 
ical probes of the upcoming decade. We adopt a Fisher 
matrix approach to forecast results achievable with a 
next generation galaxy redshift and weak lensing sur- 
vey. There are many large weak lensing and galaxy red- 
shift surveys proposed for the coming years; we choose to 



1 dx dx 



(38) 



The Fisher matrix (FM) is the inverse of the covariance 
matrix (CM) of the parameters. Working with the FM is 
particularly useful since it allows forecasted constraints 
from different surveys A and B to be easily combined 
to provide combined forecasted uncertainties from both 
surveys, so that F^^^^^ = + F^ . In the cases we 
consider here, the observable is generally a power spec- 
trum (multipole spectrum for the CMB, galaxy power 
spectrum for the galaxy redshift survey, and convergence 
power spectrum for weak lensing). 



A. Forecasting For Galaxy Redshift Surveys 



Tegmark |128j showed that for a galaxy redshift survey, 
the Fisher matrix is given by: 



(fk / d\ii P{k)\ / d\ii P(k) 



Kff(fe) 

(39) 
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The 'effective volume', "\4ff, essentially provides a 
weighting for the uncertainty in each P(/c), and (assum- 
ing constant galaxy density) is given by: 



B. Forecasting For Galaxy Weak Lensing 



nPjk) 
1 + nP{k) 



(40) 



where Vo is the volume of the survey and n is the galaxy 
density. The power spectrum that is measured is that of 
the galaxies, not the underlying dark matter. The two 
can be related by including the effects of bias and red- 
shift space distortions. The galaxy bias b is defined as 
Sg=bSm^ where Sm is the matter density fluctuation, and 
Sg is galaxy density fluctuation. Furthermore, the true 
radial position of the galaxies as measured by the redshift 
is distorted by the peculiar velocities of the galaxies. On 
large scales, coherent infall into over densities causes an 
apparent compression along the line of sight, which de- 
pends on the growth rate, fg [l29j. On small scales, large 
peculiar velocities due to thermal motions in virialised 
clusters can elongate the positions of galaxies along the 
line of sight. The biased, redshift space galaxy power 
spectrum we used is based on the model used by [130J, 
and is given by : 



P(fc,M) = 52 



1 



1 + (fc/XCr„)2/2 



1 + P{k) (41) 



where /i is the cosine of k modes parallel and perpen- 
dicular to the line of sight. The term accounts for the 
linear galaxy bias, and the ay term introduces small scale 
velocity dispersion. The / term models the large scale 
growth. The Fisher matrix can now be written in terms 
of fi and the magnitude of k. 



-1 Amin 



27Tdkdjii f 9 In P(/c,/i) 



' d \ii P{k,ii) 



VM (42) 



We choose a model to be similar to the Euclid galaxy 
redshift survey, consisting of 18,000 square degrees be- 
twee n red shift 0.5 and 2.1, in 15 independent redshift 
bins |l3l( . Due to uncertainty at large scales, the Fisher 
matrix is insensitive to the choice of /cmin- Follow- 
ing |l32( . we set /cmax so that the variance is given by 
cr'^iji^) = 0-25, and also that A^^ax < 0.2/i-^Mpc at 
higher redshifts. We assume a constant galaxy density 
of 1 X 10~^/i~^Mpc^ and a constant linear bias of 1.7. 
The volume between each bin for a survey siz e of 18,000 
square degrees was ca lcula ted according to jl33(, with 
code implemented by |l34j . Full details of the redshift 
bins is included in Appendix [A] To calculate the total 
Fisher matrix for each survey, individual Fisher matrices 
were calculated for each redshift bin, and then summed. 

For each fiducial model, a Fisher matrix was calculated 
for each of the 15 independent redshift bins in the survey 
model. These matrices were summed to obtain the Fisher 
matrix for the whole redshift survey. 



As shown in [135|, the Fisher matrix for the conver- 
gence power spectrum Pf^ is given by: 



WL 



y- 



where fsky is the fraction of the sky covered by a survey, n 
is the average number density of galaxies, and (7^>^t)^^^ 



IS 



the galaxy intrinsic rms shear. We took (7^^^)^/^ = 0.22 
(following [22]) to match the Euclid experiment, fsky = 
0.5 and n 4.18 x lO^sr"^, and the galaxy distribution 
function given by Eqn. [3ll so that n{z)dz = 1, with 
zo = 0.9. 

However, since we are interested in Euclid-like surveys, 
we used weak lensing tomography which has been proven 
to improve significantly the constraints on cosmological 
parameters (ll4j . We used 5 equal-galaxy- number bins, 
for z = — 3. To account for the photometric redshift 
uncertainties, we followed |l36| . In the i^^ bin, limited 
by Zi and ^i+i, we used the distribution ni{z) defined by 



n.(.) = in(.) (erf(^ 



erf ( 



(44) 



with = 0.03(1 + z). 

In the convergence power spectrum, this will give a 
factor Wi{x)Wj{x) instead ofW'^{x)i where Wi{x) comes 
from ^i(x), itself defined with ni{z). Instead of a single 
power spectrum Pp^ we now have a matrix Pf , where 
the (i, j) component comes from the Wi{x)Wj{x) factor. 

Equation (|43]) can be generalized by 

^;r= E(^+V2)/...)Tv(cr^^cr^g^) (45) 



where is given by = P]-^ + {Im) l^i^ij^ 
the average number of galaxies in the i^^ bin. 

This expression does not take into account various 
other systematics. For instance, intrinsic alignments of 
galaxies can cause a false positive of the lensing effect: 
the ellipticities of nearby galaxies can b e correlated be- 
cause of gravitationnal interaction |l37( . We did not in- 
vestigate other sources of errors, such as the finite size 
effects of the point spread function or its anisotropy [138J . 

The details of our weak lensing survey are given in 
Appendix [Al 



C. CMB Forecasts 

We use "FisherCodes" by Sudeep Das [139] to forecast 
results that could be obtained from the CMB. We set the 
survey specifications for a Planck-style mission based on 
[Tioj . TT, TE and EE spectra were used, to Irnax = 2000, 
with fsky = 0.8. Three channels at 100, 143, 217 GHz 
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FIG. 10: One standard deviation Fisher ellipses for rriu^ ^eff 
and fax- The full combined set of parameters also includes w, 
Qbh^ , ^ch^ , Us and As, which have been marginalized over in 
this plot. Here the weak lensing forecast (blue dashed line) 
provides the best constraints. These forecasts are for an axion 
mass of iTia — 10~^° eV. 



were used, with beam full-width half- maximum set to 
9.5', 7.1' and 5.0'. The temperature noise per pixel was 
set to 2.5, 2.2 and 4.8, and the polarization noise per 
pixel was set to 4, 4.2 and 9.8 xlO^ At /Tcmb- 



D. Results: CMB, Galaxy Redshift Survey and 
Weak Lensing Forecasts 

We first consider a cosmology with massive neutrinos 
and a fraction of CDM in axions. We assume a fiducial 
model with parameters w = —1, ft^h'^ = 0.02258, ftch'^ = 
0.1109, Us = 0.963, As = 2.3 x 10"^, = 0.055eV, 
Neff = 3.04 and fax = 0.01. Following [l40j, we also in- 
clude the optical depth to reinozation r = 0.166 and he- 
lium fraction Yhc = 0.24 as free parameters in the CMB 
FM. Similarly, galaxy bias 6=1.7 and non-linear velocity 
dispersion = 350 kms~^ were included as free parame- 
ters in the GRS FM. These parameters were marginalized 
over before the matrices were combined. CMB, GRS and 
WL Fisher matrices were calculated for this fiducial cos- 
mology for rua = 10-2^, 10-^°, 10"^^ and 10"^^ eV. As 
an example. Fisher ellipses for the rria = 10~^^ eV case 
are plotted in Fig. [10] for axion and neutrino parame- 
ters. We see the expected negative correlation between 
fax and m^y, but a positive correlation between fax and 

Neff- 

With the different sets of forecasts over a range of ax- 
ion mass, we can compare how the fully marginalised 
uncertainty on our parameters varies with axion mass. 
We can see in Fig. [11] that the uncertainty in fax from 
galaxy redshift surveys does not change appreciably with 
axion mass, as expected for the masses under consider- 



FIG. 11: Marginalized uncertainty in fax for our three ob- 
servables, evaluated for four different fiducial axion masses, 
for the cosmology ACDM+/acc + rriu. The uncertainty from 
the GRS and WL surveys does not change appreciably across 
the range of axion mass, whereas the uncertainty from the 
CMB survey decreases for higher axion masses. We cannot 
be sure if the small decrement in uncertainty in the WL survey 
at TTia = 10"'^'^ is significant. 



ation, and providing a good test that uncertainties in 
the numerical background evolution, which would vary 
with mass, do not affect the constraints. However, we 
do see some unexpected variation in the weak lensing 
constraints with mass, of 0{1) of the total uncertainty. 
We believe this is due to numerical uncertainty in the 
complex calculation for lensing with an oscillating back- 
ground, and so our results here can only be considered 
reliable to within this extra uncertainty. 

The CMB becomes more sensitive to fax for higher 
mass: this we believe to be physical. It is caused by the 
increased ISW effect at redshifts close to the surface of 
last scattering, as expected. The trend to slightly better 
constraints at higher mass has been tested and seen to be 
stable: stability requires calculations to be made at ex- 
tremely high accuracy. The required level of accuracy is 
computationally very expensive, and ensuring its stabil- 
ity has been the main limitation on models we have been 
able to forecast for. The constraining power of the CMB 
alone allows fax ~ 0.1 to be detected, which is consistent 
with the results of [85]. 

The improvement in the CMB measurement of fax is 
shown again in Fig. [T2]for the range of fiducial m^, where 
we show the Fisher ellipse with l^c^^ -These results show 
the expected positive correlation between fax and Qch'^ 
caused by the effects on equality already discussed. The 
increase in CMB constraining power at high masses does 
not reach a sufficient level to break this degeneracy. 

In Fig. [13] we show the effect of having redshift in- 
formation on the constraints from a galaxy survey. Al- 
though a measurement of the power spectrum at a single 
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FIG. 12: Fisher ellipses for our forecast CMB survey for pa- 
rameters Qch^ and fax for a range of axion mass. Note that 
if we fix fax, we obtain the same uncertainty on Qch^ for all 
the different masses. 



FIG. 14: Constraints on ris and fax, for rua = 10"*^° eV, 
showing the expected positive correlation from galaxy redshift 
surveys. At this level of constraint, no correlation is visible 
in the CMB. 
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FIG. 13: Forecasts on Qch^ and fax, for the galaxy redshift 
survey. With one redshift bin of the survey (dashed line), 
the survey cannot discriminate between a cosmology with or 
without axions, and can only rule out fax ^ 10%. When 
constraints from all 15 redshift bins are combined, detecting 
fax >0 is just possible to la. The constraints are the same 
for all axion masses. 
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FIG. 15: Constraints on Qbh^ and fax, for rUa = 10~^° eV. 
We note a negative correlation from galaxy redshift surveys, 
whereas at this level of constraint, no correlation is visible in 
the CMB.. 



redshift cannot distinguish fax = 0.01 from fax = at 
Icr, this becomes possible when all z bins are combined. 
We reiterate that this constraint is independent of axion 
mass. 

Fig.[T4lshows the expected positive correlation between 
constraints on fax and constraints on Ug from GRS. The 
constraints on Qbh'^ and fax are shown in Fig. [T5l show 
a negative correlation in GRS, and small positive corre- 
lation in WL. Both figures show no visible correlation in 
the CMB on the scale of the plot. The strong constrain- 



ing power of the CMB on ftijh'^ and relative to fax 
breaks the expected correlations at this level. 

For comparison, we also considered a cosmology as be- 
fore, but with massless neutrinos. The constraints on fax 
from all observables, along with the mass dependence, 
was approximately the same as without neutrinos. 

We finally consider an axion-free cosmology, with mas- 
sive neutrinos. The Fisher ellipse for N^ff and rrijy is 
shown in Fig. [16] for our different observables. 
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FIG. 16: Constraints on A^e// and vfii; for the axion-free cos- 
mology. 



V. DISCUSSION AND CONCLUSIONS 

Ultra-light axions and sub eV mass neutrinos are on 
the face of it very similar, when viewed in their effects on 
the matter power spectrum, and we seem hard pressed 
to distinguish them from ACDM, never mind one an- 
other, within the limits of current or future observations 
of galaxy clustering alone. There are in addition many 
other cosmological ingredients that can mimic various 
parts of such a signal and display degeneracies with ax- 
ions and neutrinos. Despite this, we have seen that in 
fact quite tight bounds could be placed on the axiverse 
by combining many precision cosmological observations, 
in particular weak lensing tomography. 

All the effects of an ultra-light species stem from two 
main sources: the effective removal of some dark matter 
from the background expansion while the field is frozen at 
redshifts greater than Zqsc^ and the suppression of struc- 
ture due to the quantum pressure of the field below the 
scale /cm, an effect very similar to neutrino free stream- 
ing. This is because the axions are having purely gravi- 
tational effects. The effects of quantum pressure on the 
matter power spectrum, growth rate, and galaxy weak 
lensing can all be qualitatively understood in terms of 
the simple step picture provided by the fitting of Tax- 
The effects on the CMB, and additional effects in the 
BAO can be understood simply from the transition in 
scalar field behaviour at Zqsc 

An important positive correlation exists between the 
fraction in axions, fax^ and the total content in CDM, 
r^c^^, since these ingredients have opposite effects in 
moving cosmological scales, such as equality and Silk 
damping. However, fax can still be quite tightly con- 
strained and this degeneracy broken. In the case of GRS 
observables alone, a fraction in axions of 1% can be dis- 
tinguished from zero, independent of mass. Redshift in- 



formation breaks the degeneracy with CDM by resolving 
growth rate effects that compliment those obtained from 
the power spectrum alone. Axions cause a distinctive 
scale dependent suppression of the growth rate at the 
scale "i 

in a manner which can be understood in or- 
der of magnitude as being due to the step in the power 
spectrum. Tax- also controls the magnitude of the 

growth rate, so that measurements of growth and clus- 
tering together can constrain CDM and axions indepen- 
dently. 

The mass independence of axion constraints from 
galaxy redshift surveys and weak lensing can be im- 
portant for the axiverse scenario, since a measurement 
will consequently constrain the existence of any and all 
species in the mass range we have investigated, where 
axions begin oscillating in the matter dominated era, 
rria ^ 10~^^ eV. Such species can be resolved, at la, 
in quantities of around 1% of the total DM content. The 
constraint is slightly worse when massive neutrinos are 
not present. The small dependence on neutrinos can be 
understood since neutrino mass adds more radiation at 
early times, and consequently moves the relevant scales, 
in this case making the axion effects slightly harder to 
resolve in the surveys considered. 

The correlation of axions with baryon content is more 
subtle, in that it is observable dependent. We expected 
a positive correlation in CMB observables: axions shift 
equality one way, baryons adding more matter, shift it 
the other, but the vastly different constraining power of 
the CMB between these observables did not make this 
visible. However, there is a negative correlation in GRS, 
where baryons have the competing effect of also sup- 
pressing structure. Similarly, we see a positive correla- 
tion between axions and the number of effective neutrino 
species. One normally thinks of extra relativistic energy 
density causing suppression of power in the CMB, but 
this is only true if all other parameters are altered so 
as to preserve equality and the angular horizon size. A 
simple increase in N^ff actually boosts CMB power, and 
competes in an opposite manner to axions. If equality 
and horizon size are fixed in the presence of axions, the 
CMB is left unchanged: axions effect the CMB purely 
through the expansion rate, the change to which is only 
slight near Zosc) making the unlensed CMB a poor ob- 
servable to look for axions. 

Constraints on axions from the CMB are subtly depen- 
dent on axion mass. Although the axions effect the CMB 
purely gravitationally via the expansion rate, they do not 
simply "switch on" to DM behaviour out of nowhere, but 
contribute as a DE term prior to this. If the DE contri- 
bution is significant near to the surface of last scattering, 
then the CMB becomes more sensitive to these axions. 
This is true for the heaviest axions we have considered, 
since their masses cause slow roll to occur near to last 
scattering, and the axions transition through this while 
contributing of order fax to the energy density. The con- 
straints on this axion mass range coming from the CMB 
are thus related to well known CMB constraining power 
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on EDE. This result is also consistent with the results 
of [SSJ, where it is this mass range that is most strongly 
constrained by the CMB, with a sharp rise in uncertainty 
outside of this range when axions are either pure CDM or 
pure cosmological constant as far as the background ex- 
pansion post recombination, and hence the CMB is con- 
cerned. However, axions are distinct from EDE in that 
they exit slow roll shortly after this, whilst EDE does 
not. An investigation into correlations between such ax- 
ions and popular models of EDE would be interesting to 
pursue. 

In all cases, we have seen the tightest constraints com- 
ing from weak lensing, and just like GRS the strength 
of this constraint depends on the photometric redshift 
measurement i.e. on tomography. Lensing tomography 
allows another measurement of the growth rate, and the 
redshift evolution of the axion suppression of small scale 
convergence power can be resolved. 

Our major result has been to show that with current 
and next generation galaxy surveys alone it should be 
possible to unambiguously detect a fraction of dark mat- 
ter in axions of the order of a few percent of the total. 

In conclusion, we have seen that even such a simple 
ingredient as a light, non-interacting scalar field, such as 
the axion, can lead to interesting cosmology. The effects 
are sometimes subtle, and minor, making such an ingre- 
dient hard to spot. When these constraints are indepen- 
dent of the mass of the scalar field, they will be important 
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Appendix A: Forecast Survey Details 

In this appendix we include details of the surveys. The 
central redshift, volume, and /cmax for each of the 15 red- 
shift bins in our galaxy redshift survey are given in Ta- 
ble [III 

For the lensing survey, we define the redshift bins in 
Table lEI) and give the survey parameters in Table HVl 
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